Counting the number of correlated pairs in a nucleus 
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Abstract 

We suggest that the number of correlated nucleon pairs in an arbitrary nucleus can be estimated 
by counting the number of proton-neutron, proton-proton, and neutron-neutron pairs residing in a 
relative S state. We present numerical calculations of those amounts for the nuclei 4 He, 9 Be, 12 C, 
27 Al, 40 Ca, 48 Ca, 56 Fe, 63 Cu, 108 Ag, and 197 Au. The results are used to predict the values of the 
ratios of the per-nucleon electron-nucleus inelastic scattering cross section to the deuteron in the 
kinematic regime where correlations dominate. 

PACS numbers: 25.30.Rw,25.40.Ep,24.10.Jv,24.10.-i 



The nucleus is a prototype of a dense quantum liquid with a high packing fraction [T]. 
Naively one could expect some severe medium effects for the nucleons. Several experimental 
investigations confirmed the robustness of the nucleons. This is for example reflected in 
the successful use of the Impulse Approximation (IA) in nuclear reaction theory. In the 
IA the bound and free nucleon properties (charges, magnetic moments, form factors) are 
considered identical. A few experiments, however, found indications for medium-modified 
nuclear properties. Recent 4 He(e, e'p) measurements [2], for example, could be described 
after implementing medium-modified proton form factors. Also in comparing deep inelastic 
scattering cross sections with those on the deuteron, one finds that under some kinematics 
conditions the naive scaling ratios do not hold. This observation is known as the EMC 
(European Muon Collaboration) effect [3] and indicates that under selected kinematics the 
whole of the nucleus appears to be more than the sum of its constituents. 

Recently, it was suggested [I] that the magnitude of the EMC effect can be predicted 
from the knowledge of the measured a 2 (A/D) coefficients. The a 2 (A/D) coefficients are 
defined as 



where a A (xb, Q 2 ) is the inclusive (e, e') cross section for the target nucleus A at a particular 
four-momentum transfer Q 2 and Bjorken 1.4 < Xb = < 2 (M is the nucleon mass, 
and uj the energy transfer). The observed plateau in the measured xb dependence of a 2 for 
1.4 < xb < 2 is a strong indication for scattering from a correlated nucleon pair [5j E]. As a 
matter of fact, the a 2 coefficients can be interpreted as a measure for the effect of short-range 
correlations (SRC) in the target nucleus A relative to deuteron D. In this paper, we suggest 
a technique that allows one to estimate the number of nucleon pairs prone to SRC in an 
arbitrary nucleus A(N, Z). We use these estimates to predict the values of the coefficients 



A time-honored method to quantify the effect of correlations in classical and quantum 
systems is the use of correlation functions. The latter encode those portions of the system 
that depart from mean-field behavior. The realistic (correlated) wave functions | are con- 
structed by applying a many-body correlation operator to the mean-field Slater determinant 
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The Q reflects the full central, spin and isospin dependence of the nucleon-nucleon force but 
is dominated by the central and tensor correlations 
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where g c (7" 12 ) an d ft T ( r i2) ar e the central and tensor correlation function, S12 the tensor op- 
erator and S the symmetrization operator. The correlation functions g c and ft T determine 
the radial dependence and magnitude of the correlations. Over the last couple of decades, 
various many-body calculations adopting a plethora of techniques [7HTU] have made predic- 
tions for the correlation functions g c and f tT . These calculations, confirmed the following 
robust features. First, the two-nucleon correlations represent a local property. This implies 
that the correlations are universal or only weakly A dependent [UJ. This means that g c 
and ftr are very much confined to the bulk part of the nuclear density and only depend 
on the inter-nucleon distance. The universality property implies that the f tr (r^) Sijfi.fj 
correlation operator in a nucleus A is not very different from the one that mixes the 3 Di and 
3 Si wave-function components in deuterium [Sj. Second, it was observed that for moderate 
relative pair momenta (300 < k 12 < 600 MeV), the effect of the tensor correlations is dom- 
inant [T2l [T3] . As the Su exclusively affects nucleon pairs in a spin S — 1 state, it makes 
the proton-neutron (pn) correlations to dominate at moderate values of the relative pair 
momentum. We stress that the universality property does not imply that the correlation 
functions g c and f tT are insensitive to model assumptions. The correlation functions depend 
on the choice of the Hamiltonian, for example. Indeed, a softer Hamiltonian (implying less 
correlated wave functions) will require other correlation functions than a hard Hamiltonian 

Upon computing the response of the nucleus to some one-body operator Q = J2t=i W> 
into lowest order the effect of the correlations can be implemented by means of an effective 
transition operator which includes the effect of the correlations [El [12] 
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FIG. 1: The distribution of the relative quantum numbers I = S,P,D,F,G,H,I,> J for (a) 
the proton-neutron pairs, (b) the proton-proton pairs, and (c) the neutron-neutron pairs for the 
various target nuclei. For the proton-neutron pairs there are contributions from 1 Sq(T = 1) and 
3 Si(T = 0). The contribution from the 1 So(T = 1) is indicated by the dashed line. Results are 
obtained in HO basis with tka(MeV) = 45,4" 5 - 25A~I. 



Obviously, through the correlations a typical one-body operator (like the 7* - nucleus in- 
teraction in the IA approximation) receives two-nucleon contributions which are completely 
determined by the product of the correlation functions and the one-body operator. The 
nucleon-nucleon correlations are very local and will only affect nucleon pairs which are 
"close". Accordingly, the correlation operators — g c (r^) +t act as projection operators 
and will almost exclusively affect nucleon pairs that reside in a relative S state. 

We suggest that the significance of two-nucleons correlations in a certain nucleus A(N, Z) 
is proportional to the number of relative S states. In order to compute this number, a 
coordinate transformation from (r*i, r 2 ) to (?\2 = fi — r 2 , R = ri + r2 ^ is required. The single- 
particle states in the Slater determinant | \P ) are denoted by a a = (n a l a j a m a t a ), where 
t a — ±| is the isospin quantum number. In a harmonic-oscillator (HO) basis the normalized 
and antisymmetrized two-nucleon wave functions can be written as 
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where T (S) is the total isospin (spin) of the pair. Further, \nl) (|JVA)) is the relative (center 
of mass, cm.) pair wave function. The explicit expression for the coefficient C can be found 
in Eq. (20) of Ref. [T5]. With the aid of the above expression ^ one can project two-nucleon 
states in (ri,^) on nucleon states in (?\2, Rj and determine for each pair (a a otb) of shell- 
model states the weight of the various relative (nl) and cm. (NA) quantum numbers. For 
two-nucleon states in a non-HO basis, one can obtain the weights of the various (nlNA) 
combinations by expanding the single-particle wave functions in a HO basis. 

We have computed the C coefficients for all target nuclei A either for which the a-i(A/D) 
coefficient has been published or for which one may expect data in the foreseeable future 
|16j . The Slater determinant is constructed by filling the single-particle states as they are 
determined in the nuclear shell model. We denote the Fermi level for the proton and neutron 
single-particle states as a p F and a F . The quantity 

^2 ("a" 6 ' J R M R WaOtb] JrM r ) = NZ , (6) 

JrMr a a <a p F a b <a™ 

determines exactly the number of proton-neutron pairs. Similar expressions hold for the 
number of proton-proton ( ^^f^- j and neutron-neutron pairs ^ N{N-i) j After inserting 
the right-hand side of Eq. (|5| in the above expression, one can compute how much of each 
combination \(nl, NA)LMl, SMs,TMt) of pair quantum numbers contributes to the total 
number of pairs. Here, we are particularly interested in the quantum numbers (nl) of the 
relative wave function. We denote the relative orbital angular momentum I = 0,1,2,... 
as S, P, D, . . .. The numerical calculations get increasingly more time consuming as A in- 
creases due to the combinatorics of all possible shell-model pairs. The accuracy of the 
numerical calculations can be checked against the normalization condition of Eq. ([6]). In 
Fig. [T] we display the relative contribution of the various I to the pair wave functions 
\(nl,NA)LM L , (f§) SM S , (§§) TM T ) for the nuclei 4 He, 9 Be, 12 C, 27 A1, 40 Ca, 48 Ca, 56 Fe, 
63 Cu, 108 Ag, and 197 Au. It is obvious that with increasing A a smaller fraction of the nucleon 
pairs resides in a relative S state. Whereas, for 12 C about 50% of the pn pairs has I = for 
the heaviest nucleus 197 Au this is a mere 10%. Accordingly, with increasing A, a smaller and 
smaller fraction of the nucleon-nucleon pairs will be prone to correlation effects. In addition, 
there is a strong isospin dependence as the fraction of the proton-neutron pairs residing in 
a relative S state is substantially larger than for proton-proton and neutron-neutron pairs. 
Naively, one could expect that the number of correlated pn (pp) pairs in a nucleus scales 
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FIG. 2: (Color online) The computed number of pp, nn and pn pairs with / = 0. For pn we 
discriminate between 3 Si(T = 0) and 1 So(T = 1). Unless indicated otherwise the results are for 
a HO basis. For the 3 Si(T = 0) pn pairs also the predictions in a WS basis are shown. The 
parametrizations for the WS potentials are from Ref. |17j . 

like NZ (^p 11 ) ~ A 2 . As illustrated in Fig. ffl our calculations rather indicate that the 
number of pairs that are prone to correlation effects follows a power law ~ ^i-44±o.oi a 
matter of fact, we find that the power law is very robust. Calculations with Woods-Saxon 
(WS) wave functions, for example, result in a computed number of S states that is very 
close (order of one percent) to the HO predictions. The N — Z asymmetry is reflected in an 
unequal number of pp, nn, and pn 1 Sq(T = 1) pairs. We stress that the ratio of the nn to pp 
1 Sq(T = 1) pairs can be considerably smaller than predicted by naive ^(z-i) comDma t° r i cs - 
For Au, for example, one expects a ratio of 2.24 whereas the data of Fig. [2] lead to 1.77. 

Now, we wish to connect the number of pairs with I = with the measured values of 
a,i(A/D). In an inclusive A(e,e') process the correlated part of the electron-nucleus (eA) 
response (corresponding with the last two terms in Eq. Q) can be probed by selecting 
events 1.4 < xb < 2. The magnitude of the response is proportional with a product of 
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FIG. 3: (Color online) The computed values for the a%(A/D) for various nuclei. The data are from 
Refs. [5J (SLAC), [6j (JLAB Hall B) and [18J (JLAB Hall C). The triangles denote the theoretical 
predictions obtained with the Eq. ([7]). 

two terms. First, the number of pairs that are prone to SRC, and, second, the value of 
the correlation functions evaluated at the relative momentum of the pair. Indeed, as is 
pointed out in Refs. [T^J ED] in the kinematical regime where correlations are probed, the 
eA response obeys ~ F{P)<j e ^N{ki 2 ) 1 where P is the cm. momentum of the correlated 
pair on which the absorption takes place and F(P) is the corresponding cm. distribution 
(the combination F(P)a e NN{ki2) is referred to as the decay function in Ref. [IS]). The 
a e NN stands for the elementary cross section for electron scattering from a correlated NN 
pair. The u e NN contains the Fourier-transformed correlation functions gdjkiz) and ft T {ki2) 
evaluated at the relative momentum fc 12 of the pair. An analytic expression for a epp can be 
found in Ref. [20]. It is worth stressing that given the kinematics, there are two possible 
values of corresponding with photoabsorption on nucleon "1" and photoabsorption on 
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Ref. [5] Ref. [6] Ref. 0] Ref. [S] Eq. (7) 



4 He 
9 Be 

12 C 

27 Al 
40 Ca 
48 Ca 

56 Fe 

63 Cu 
108 Ag 

197 Au 



3.3 ±0.5 3.80 ±0.34 3.60 ±0.10 2.4 

4.08 ±0.60 3.91 ±0.12 2.8 

5.0 ± 0.5 4.75 ± 0.41 4.75 ± 0.16 3.3 

5.3 ±0.6 5.13 ±0.55 4.4 

5.44 ±0.70 5.2 

5.2 

5.2 ±0.9 5.58 ±0.45 5.4 

5.21 ±0.20 5.6 

7.29 ±0.83 6.3 

4.8 ±0.7 6.19 ±0.65 5.16 ±0.22 7.0 



TABLE I: The a,2(A/D) values for various nuclei. The data from direct measurements of the 
nucleus to deuteron cross sections are from Refs. [5] (SLAC), [6] (JLAB Hall B) and p2] (JLAB 
Hall C). The values of Ref. [I] are phenomenological extractions based on the measured EMC data 
and the observed linear correlation between the magnitude of the EMC effect and the measured 02 
scaling factor. The quoted values of Ref. [IB] are the raw ratios. Ref. [THJ also contains corrected 
values for 02 which are about 15% smaller. 

nucleon "2" of the pair. The dominant contribution to the inclusive A(e, e') cross section 
for 1.4 < x# < 2 stems from pairs with kp < ^12 ^$ 2&f, with kp the Fermi momentum. In 
that momentum region, the g c {ku) is substantially smaller than ftriku), which causes the 
tensor correlated pn pairs to dominate [21] [22] [23] . 

The universality of the tensor correlations, which translates to the weak A dependence 
of ft T (ki 2 ), allows one to assume that the cross section a epn for electron scattering from a 
correlated proton-neutron pair in the nucleus will almost equal the one for electron scattering 
from the deuteron, provided that the cross sections are evaluated at equal values of the high 
relative momentum kyi of the pair. In a symbolic way, this feature can be expressed throught 
the scaling relation a epn (ku) ~ cr e p,{ki2). This property is related to the fact that at high 
momenta the nuclear momentum distributions n A {k) are very much like scaled deuteron 
momentum distributions: n A {k) ~ C A n D (k), where C A is a measure for the amount of pn 
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correlations in A [TT] |24j . 

With the above-mentioned scaling relation <J epn (ki 2 ) ~ cr e D (^12) valid at high relative 
momenta, one can transform the ratio of of Eq. ([I]) (the per-nucleon electron-nucleus inelastic 
scattering cross section to the deuteron) into the form 



where the integrations extend over those parts of the phase space (PS) which are com- 
patible with 1.4 < xb < 2. The quantity B^ Q (A) is the number of pn pairs in a rel- 
ative \n, I) state with the quantum numbers of the deuteron, 3 Si(T = 0). One can es- 
timate the B^ (A) from Eq. ^ by combining the computed coefficients for all possible 
|(n = , 1 = 0, NA)LM L , S = 1M S , T = 0M T = 0) combinations. In Fig. [2] we have summed 
over all possible n to obtain the total amount of I = states. For all target nuclei, the n = 
contribution dominates, but its relative importance decreases with growing A. The n = 
represents 100% of the I = pn states for 4 He, about 80% for the medium-heavy nuclei 
(Ca, Fe, Cu), 70% for 108 Ag, and 62% for 197 Au. Pairs residing in a \n 7^ 0, 1 = 0) state have 
a much smaller chance of being "close" than their \n = 0,1 = 0) counterparts and are less 
prone to SRC effects. We assume that only |n = 0, 1 — 0) proton-neutron pairs contribute 



The cm. motion of the pair in finite nuclei (absent in the deuteron) and the imposed 
conditions in xb make that a fraction of the correlated proton- neutron pairs are not counted 
in the A(e,e') signal in the numerator of Eq. Q. The F(P) for a nucleus A can be reli- 
ably computed in a mean-field model. Indeed, the 12 C(e,e'pp) measurements of Ref. [25] 
determined F(P) over a large P range and observed it to be compatible with a mean-field 
prediction. We have performed Monte-Carlo simulations in order to determine the correc- 
tion factor j ps dPF(P) for all nuclei which are considered here. We find that for A > 4 
about 25% of the correlated pairs are excluded from the experimentally scanned phase space 
due to the cm. motion of the correlated pair. From the simulations we observed that the 
correction factor is only slightly mass number dependent. With this correction factor, the 
Eq. ([7]) allows us to make predictions for the 02 (A/D). The predictions are contained in 
Fig. [3] and Table |T] and compared with experimental data. One striking observation from our 
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calculations, is that the predicted a 2 (A/D) for 40 Ca and 48 Ca are identical and equal to 5.2. 
On the basis of naive NZ combinatorics one may have expected a 30% difference between 
the two. For heavier target nuclei, the data seem to suggest that the a 2 (A/D) coefficient 
saturates. Our calculations predict a strong linear rise in the A dependence of the a 2 for 
A < 40. At higher A one enters a second regime with a much softer linear rise with A. Our 
calculations increase linearly with log (A) and tend to underestimate the data at low A and 
overestimate the data for the heavier nuclei. Final state interactions, for example, which 
are neglected in this work, may induce some additional A dependence in the a 2 ratio |18j . 
It is clear that more data are needed to establish the situation at large A. The observed 
phenomenological linear relationship between the scaling factor 02 and the magnitude of the 
EMC effect pi] gives a 2 = 7.29 ± 0.83 for Ag and a 2 = 6.19 ± 0.65 for Au, values that are 
not inconsistent with our results. 

In conclusion, we suggest that the number of correlated pairs in a nucleus is proportional 
with the number of relative S two-nucleon states. We find this number to obey a power 
law dA 1A4±om with d = 0.39 ± 0.02 (d = 0.13 ± 0.01) for T = (T = 1) proton-neutron 
pairs. The power law is robust in that it is independent of the choices made with regard to 
the single-particle wave functions. We have used the computed amount of T = pn pairs 
to predict the value of the measured a 2 (A/D) coefficients, which provide a measure of the 
number of correlated pairs in the target nucleus A relative to the deuteron. The observed 
power law in the number of relative S states translates to a linear increase of a 2 with log(A). 
We observe that our predictions are not inconsistent with the trend and magnitude of the 
data, lending support to our suggestion. 

This work was supported by the Fund for Scientific Research Flanders. We are grateful 
to Kris Heyde for useful discussions and comments. 
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